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ABSTRACT 

Accretion flows around black holes generally result in mass-outflows that exhibit irreg¬ 
ular behavior quite often. Using 2D time-dependent hydrodynamical calculations, we 
show that the mass-outflow is unstable in the cases of thick accretion flows such as the 
low angular momentum accretion flow and the advection-dominated accretion flow. 
For the low angular momentum flow, the inward accreting matter on the equatorial 
plane interacts with the outflowing gas along the rotational axis and the centrifugally 
supported oblique shock is formed at the interface of both the flows, when the vis¬ 
cosity parameter a is as small as a ^ 10 -3 . The hot and rarefied blobs, which result 
in the eruptive mass-outflow, are generated in the inner shocked region and grow up 
toward the outer boundary. The advection-dominated accretion flow attains finally in 
the form of a torus disc with the inner edge of the disc at 3 R g < r < 6 R g and the 
center at 6 R g ^ r ^ 10i? g , and a series of hot blobs is intermittently formed near the 
inner edge of the torus and grows up along the outer surface of the torus. As a result, 
the luminosity and the mass-outflow rate are modulated irregularly where the lumi¬ 
nosity is enhanced by 10 — 40% and the mass-outflow rate is increaed by a factor of few 
up to ten. We interpret the unstable nature of the outflow to be due to the Kelvin- 
Helmholtz instability, examining the Richardson number for the Kelvin-Helmholtz 
criterion in the inner region of the flow. We propose that the flare phenomena of Sgr 
A* may be induced by the unstable mass-outflow as is found in this work. 

Key words: accretion, accretion discs - black hole physics - hydrodynamics - Galaxy: 
centre. 


1 INTRODUCTION 

The observations as well as numerical simulations around 
black holes demosntrate regular, irregular and/or chaotic 
time variation of emergent radiations and mass-outflows 
over a wide range of time scales. The origins of such 
variable emissions and mass-outflows from the vicinity of 
the black holes are also numerous. The accretion luminos¬ 
ity depends on the adopted accretion models and its ac¬ 
cretion rates. When the accretion rate is not too high, 
the luminosity is directly proportional to the accretion 
rate and can be described successfully by the standard 
thin-disc model (Shakura -Sunyaev, hereafter S-S, model) 
dShakura fe Sunvaevll 19731) . The standard model have been 
proven to be successful in studying the cataclysmic vari¬ 
ables and later for neutron stars and black holes. However, 
for highly luminous accretion discs whose luminosity exceeds 
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the Eddington luminosity as is found in the black hole can¬ 
didate SS 433 or for very low luminous discs compared with 
the predicted accretion rate as the case of the supermassive 
black hole candidate Sgr A*, the S-S model is not valid and 
alternative thick disc mo dels, such as the slim disc model 
ll Abramowicz et alj 198 81). the advection-d ominated accre¬ 
tion disc model ( Naravan fe YilPl 9941. Ii9 95l) and the low an¬ 
gular momentum disc model (Moscibrodzka. Das fe Czernvl 


I2OO6I : ICzernv fe Moscibrodzka 120081 ). have been developed. 


Numerous studies of the hot thick accretion flow around 
the black holes have been carried out both theoretically and 
numerically to promote the insights of the black hole physics. 
One of the important results of the hot accretion flows is that 
the inflow and outflow rates in the accretion flow follow a 
power-law function of radius. To explain the inflow-outflow 
rate versus radius relation, tw o types of models of the adia¬ 
batic inflow-outflow so lution dBlandford fe Begelmanl|l999l . 
1 20041 : I Begelmaril 120121) and the convection-dominated ac¬ 
cretion flow model dNaravan, Igumenshchev fe Abramowicz! 
l2000l : iQuataert fe Gruzinov 2000l : Abramowicz et al.l 20021 : 
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lleumenshchevlliood 'l have been proposed. Furthermore, the 
recent hydrodynamical (HD) and magneto-hydrodynamical 
(MHD) simulations result in new findings of such 
outflows, winds and jets in the field of hot accre- 


tion flows around the black holes (Narayan et al. 2012!: 

Yuan. Bu & Wu 201 4 Yuan. Wu & Bui 2012 

: Bu et al. 

20lJ Li. Ostriker & Sunvaevll2013l: Yuan et al. 

201.5l). The 


details o f the g e omet ricall y thick hot accretion flow s are re¬ 
ferred to lYuanl ll201ll ) and I Yuan fc Naravanl ll20l4) . 


The hot accretion flows show unstable behaviors of the 
outflows in most cases of the time-dependent numerical sim¬ 
ulations. We have examined the thick accretion flows with 
low angular momentum and found that the outflows from 
the outermost boundary are always unstable, accompanying 
irregularly er uptive mass ejections llOkuda fe Moltenill2012l : 
IOkudall2014l l. Although we interpreted the unstable mass- 
outflow to be induced by the irregular radial oscillations of 
the centrifugally supported shock, the intrinsic origin of the 
unstable nature of mass-outflows is still an open question. In 
this paper, focusing on the outflows from the hot accretion 
flow, we examine the unstable nature of the outflows in terms 
of the Kelvin-Helmholtz instability, using time-dependent 
2D hydrodynamical simulations. Finally, following the re¬ 
sults of the unstable mass-outflow, we suggest a possible 
model which may explain the flaring events of Sgr A*. 


2 NUMERICAL METHODS 
2.1 Basic equations 

We consider a supermassive black hole with mass M = 
4 x 10 6 Mq which is estimated for the supermassive black 
hole candidate Sgr A*. We also assume a two-temperature 
plasma with ions being much hotter than electrons. The ac¬ 
creting gas is assumed to be very rarefied and optically thin. 
The set of relevant time-dependent equations is given in the 
spherical polar coordinates (r,£,(/?): 

^ + div(pu) = 0 , ( 1 ) 
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where v = (v,w,v v ) are the three velocity components, p 
is the density, £i and £ e are the specific internal energies 
of the ion and electron, pi and p e are the gas pressures 
of the ion and electron, S = (S r ,S^,S v ) denotes the vis¬ 
cous stress tensor and <f> = (S v) v is the viscous dissi¬ 
pation rate per unit mass. The full expression of the vis¬ 
cous stress tensor S and the dissipation rate $ are given in 


lOkuda, Fuiita fe Sakashital dl997l l. Here, the kinematic vis¬ 
cosity v is given by 


v — ac a min [H p , H], 


(7) 


where c s , H p and H are the local sound speed, the pressure 
scale height and the disc thickness, respectively, and a is 
a dimensionless viscosity parameter which is usually taken 
to be in the range 1(D 3 ^ a < 1. The pseudo-Newtonian 
potential is adopted in the momentum equation. q el , g br 
and Qsyn are the energy transfer rate from ions to elec¬ 
trons by Coulomb collisions, the cooling rate of electrons 
by electron-ion and electron-electron bre msstrahlungs and 
the synchrotron coo l ing rate, respectively dEsin et al.lll996l : 
iNaravan &; Yil 1 19951 : IStepnev fe GuilberH 1 19831 ) . In calcula- 
tion of q sy n, we give the magnetic field B at each radius, 
assuming that the ratio /3 of the magnetic energy density to 
the thermal energy density is constant throughout the re¬ 
gion. We neglect the radiation transport assuming that the 
flow is optically thin. 

In the equations (5) and ( 6 ), we assume that all of the 
viscous dissipation energy is given to ions because an ion 
particle is much heavier than an electron particle and it is 
poorly known how the viscous dissipation energy is shared 
between the ions and electrons. However, as a representa¬ 
tive case, we examined what result is obtained when half 
of the dissipation energy is shared to the electron and ob¬ 
served that the total luminosity is altered maximumly by a 
small factor of 1.5 due to the larger electron temperature in 
the intermediate region on the equator, although the mass- 
outflow rate is not altered largely. This indicates that the 
above assumption is reasonable for our present study. 

Recent MHD simulations show that r — <p component 
of the viscous stress is dominated over other components. 
However, we used the true kinematic viscous stress which 
includes all components of the stress, in expectation of en¬ 
hancement of the outflow activity through the angular mo¬ 
mentum transfer and the viscous dissipation. From some 
examinations of the effects of viscous stress components, we 
find that in both cases B and C the global features of the 
total luminosity and the mass-outflow rate are not altered 
largely even if only the r — ip component of the stress is used. 

In the present two-temperature plasma with the ion 
temperature Ti and the electron temperature T e , the adi¬ 
abatic indices for ions ( 71 ) and electrons ( 7 e ) possess dif¬ 
ferent values depending on their temperature states. In re¬ 
ality, we use 7 e is 1.6 for kT e ^ m e c 2 , equivalently, T e ^ 
5.9 x 10 9 K(= T c ) where electrons are non-relativistic and 7 e 
is 4/3 for kT e ^ m e c 2 where electrons are relativistic. Here, 
T c is the critical electron temperature at which electrons in 
the non-relativistic state change into the relativistic state. 
For ions, 71 is taken to be 1.6 throughout the region b ecause 
ions remain in a non-relativistic state of kTi ^ m p c 2 (iFukud 
1 19861) . 

The set of partial differential equations (1)—( 6 ) is nu¬ 
merically solved by a finite-difference method under ad¬ 
equate initial and boundary conditions. The numerical 
schemes used are basically the same as that described by 
|Ok udal201 j l. The methods are based on an explicit-implicit 
finite difference scheme. 
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2.2 Modeling of the thick accretion flow 

We focus here on the geometrically thick and optically thin 
accretion discs. The low angular momentum flow and the 
advection-dominated accretion flow around the black holes 
belong to this category. Using the flow parameters of the 
thick accretion flow, we examine the time-dependent behav¬ 
iors of the outflow which is ejected from the outer bound¬ 
ary. We consider three representative cases of the accretion 
flow with no angular momentum (case A), the low angu¬ 
lar momentum flow (case B) and the advection-dominated 
accretion flow (case C), respectively. In the previous work, 
we have studied t he low angular momentum flow mod el in 
the inviscid limit (lOkuda fe MoltenHl2012l : lOkudalfeo Mh but 
we intend to re-examine here considering the effect of vis¬ 
cosity. The case A is studied for initial flow with flow pa¬ 
rameters same as case B. For case C, we start with the 
self-si milar solutions of the advetion-dominated accretion 
flow dNaravan fc Yi|[l995l : lYuan fc Naravan|[20l4l . Follow¬ 
ing this, the self-similar solutions are given by 


v « -0.367 ar _1/2 , (8) 

A « 0.292 r 1/2 , (9) 

c s 2 « 0.156 r"\ (10) 

p = 1.05 x 10 -4 a” 1 m _1 mr^ 3 ^ 2+s gem -3 , (11) 


where A is the specific angular momentum, s is the index pa¬ 
rameter. Here, r is measured in units of the Schwarzschild 
radius R s , v and c a are measured in speed of light c and 
A is given in units of cR g . The mass of the black hole m 
and the accretion rate m are expressed in unit of solar mass 
and the Eddington rate of 1.4 x 10 18 gm s -1 , respectively. 
In this self-similar solution, we consider a constant mass 
accretion rate with the index s = 0. Although recent HD 
simulations of the thick hot accretion flow over a wide range 
of r ~ 10 4 i? g (lYuan. Wu fe Bull2012l l reveal slightly differ¬ 
ent relations compared to the above self-similar solutions 
(i.e. p oc r -0 ' 65 and r -0 ' 85 for a =0.001 and 0.01, respec¬ 
tively), eventually it does not affect the present study of the 
mass-outflow. In Table 1, we show the flow parameters for 
cases A, B and C, where the flow variables at the outer disc 
boundary are illustrated. The electron temperatures at the 
outer boundary are approximately taken to be equal to the 
ion temperatures. 

2.3 Initial and Boundary Conditions 

The initial discs of cases B (also A) and C are given by 
ID sol utions of the low angular momentum accretion flow 
model (IOkudall2Q14l f and the self-similar solutions (8)-(ll) 
of the advection-dominated accretion flow, respectively. Fur¬ 
thermore, assuming the density distribution of p ~ r -1 , we 
construct an initial rarefied atmosphere around the discs to 
be approximately in radially hydrostatic equilibrium. 

Physical variables at the inner boundary, except for the 
velocities, are given by extrapolation of the variables near 
the boundary. However, we impose limiting conditions that 
the radial velocities are given by a fixed free-fall velocity and 
the angular velocities are zero. On the rotational axis and 
the equatorial plane, the meridional tangential velocity w is 
zero, and all scalar variables must be symmetric relative to 



Figure 1. Time variations of total luminosity Lt (erg 1 ) and 
mass-outflow rate M 0 u t in case A with the viscosity parameter a 
= 0.0 (solid lines) and 0.1 (dotted lines). 


these axes. The outer boundary at r = r ou t is divided into 
two parts — one is the disc boundary through which matter 
is continuously entering from the outer disc with a constant 
accretion rate M and the other is the outer boundary region 
above the accretion disc. We impose free-floating conditions 
on this outer boundary and allow matter to eject in the form 
of outflow where any inflow is prohibited. All flow variables 
at the outer boundary of the disc are kept constant always. 


3 THE NUMERICAL RESULTS 

3.1 The thick accretion flow with no angular 
momentum 

Case A is the simplest case of thick accretion flows. The 
initial flow variables of case A are given by ID solution of 
the low angular momentum flow model (case B) in Table 
1 but the angular momentum is neglected here. After the 
time integration, the upstream accreting gas falls towards 
the gravitational center and soon reaches close to the in¬ 
ner edge of the disc. Though most of the accreting gas is 
primarily swallowed into the black hole, a part of the gas is 
ejected upward along the z-axis due to the high temperature 
near the inner edge. Since the ejected gas interacts with the 
downward accretion gas with considerable amount of hori¬ 
zontal velocity component, an oblique shock is formed along 
the z-axis and is elongated in the region of z ~ 80R g and 
r ~ 30-Rg at t ~ 10 5 s. Eventually, the upward moving gas 
in the shocked region is escaped from the outer boundary. 
As a result, a few to several percent of the input accreting 
matter is ejected through the outermost boundary. 

Fig. 1 shows the time variations of the total lumi¬ 
nosity Lt integrated in the computational domain V and 
mass-outflow rate M 0 ut ejected from the atmospheric outer 
boundary S for case A. Here, Lt and M ou t are given by 

Lt = j (qbr + <?syn)dV, (12) 

and 
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Table 1. Flow parameters as specific angular momentum A, adiabatic index of ions 7 ;, mass accretion rate M (M@ yr —1 ), radial velocity 
fout, density p, ion temperature T;, relative disc thickness h/r, viscosity parameter a and ratio 3 of the magnetic energy density to the 
thermal energy density at Rout = 200 Rg in the two-temperature model. 


case 

A 

7i 

M (M e yr" 1 ) 

~Vout/ c 

P 

T ; (K) 

h/r 

a 

p 

A 

0.0 

1.6 

4.0 X 10“ 6 

5.06 x 10 -2 

5.61 x 10“ 19 

1.292 x 10 9 

0.43 

0 . 0 - 0.1 

10“ 3 

B 

1.35 

1.6 

4.0 x 10“ 6 

5.06 x 10 -2 

5.61 x 10 -19 

1.292 x 10 9 

0.43 

10 “ 3 -0.01 

10- 3 

C 

4.13 

1.6 

1.5 x 10“ 6 

2.60 x 10“ 2 

4.14 x 10~ 18 

2.082 x 10 9 

0.50 

0 

1 

CO 

0 

10“ 3 


Afout — J Poll t ttout d*S, (13) 

where the suffix ‘out’ shows the value at the atmospheric 
outer boundary. The solid and dotted lines represent the 
cases with a = 0.0 and 0.1, respectively. We find that the lu¬ 
minosities and the mass-outflow rates vary quite irregularly. 
In particular, the variation for first one is less than 20 % and 
for second one is by a factor of few. Here, the luminosities 
are mainly dominated by the bremsstrahlung emission. The 
mass-outflow rates increase with the increasing viscosity and 
conversely the luminosities decrease especially in the case of 
the maximum viscosity with a = 0.1. Due to the viscous 
heating, ion temperature becomes higher and therefore, the 
shock wave near the z-axis is pushed outward. Hence, the 
outflow over the oblique shock is enhanced as seen in Fig. 2. 
On the other hand, the electron temperature behind the 
shock becomes lower when the shock is pushed away. As 
the dominant energy emission process is due to the free-free 
emission which mostly depends on the electron temperature 
behind the shock, the luminosity evetually becomes lower 
by a small factor. 

We wonder why the mass-outflow rates are considerably 
modulated chaotically in spite of the small modulations of 
the luminosity. To understand the outflow activity, it is use¬ 
ful to follow the animation of the time-dependent density 
and temperature contours. Fig. 2 shows the snapshots of 
the animation over a typical cycle in the time variation of 
the mass-outflow rate, where the temperature contours at 
t/10 6 = 2.899, 2.929, 2.960, 2.990, 3.005 and 3.036 s are 
shown. The large mass-outflow activity begins at the elon¬ 
gated region of the shocked surface. Firstly, a hot and rar¬ 
efied blob with small size is formed at the upper edge (z ~ 80 
and r ~ 40) of the shocked region. The hot blob grows up in 
the outward direction and starts disintegrating into several 
small blobs. Ultimately, these blobs further continue to ex¬ 
pand outward. The maximum mass-outflow rate is attained 
at the time of Fig. 2(d). Thus, the unstable mass-outflow is 
observed even in the case of the thick accretion flow with no 
angular momentum. 


gion and consequently the luminosity and mass-outflow rate 
are modulated within a factor of few and several, respec¬ 
tively. In this work, we study the low angular momentum 
flow model including viscosity and observe that the charac¬ 
teristics of the low angular momentum flow are not altered 
largely compared to the inviscid flow, when viscosity param¬ 
eter a is small as ~ 1CP 3 . Similarly to case A, an oblique 
shock surface is formed due to the interactions between the 
upward moving gas and the downward accreting gas. Since 
the flow is rotating around the z-axis, the shock surface on 
the equatorial plane is pushed away in the outward direc¬ 
tion due to the centrifugal force and is oscillated irregularly 
at r/R g ~ 20 - 50. During shock oscillations, a small hot 
blob is occurred in the shocked region near the equatorial 
plane, grows upward, undergoes to become maximal at the 
upper edge of the shocked region and finally attains the outer 
boundary. The evaluation of this hot blob seems to be suit¬ 
able to explain the episodic mass-outflow activity. 

Fig. 3 represents the time variations of the total lumi¬ 
nosity Lt erg s -1 , the mass-outflow rate M ou t, the mass- 
inflow rate Mi n at the inner edge and the shock position 
R s /R g on the equatorial plane for case B with a = 1CU 3 . 
The mass-outflow rate attains 10 % of the input accretion 
rate and varies by a factor of 2 - 3 whereas the luminos¬ 
ity is modulated by a factor of 1.5. The shock position on 
the equatorial plane oscillates irregularly around r ~ 35 Rg. 
The large mass-outflow activity occurs roughly in regular 
intervals of the order of few days. Fig. 4 shows the contours 
of ion temperature T; at t = 1.83 x 10 6 s for case B with 
a = 1CU 3 . The time is indicated by ‘+’ on the curve M ou t 
in Fig. 3. The shock surface is illustrated by the black thick 
lines which are located at r ~ 35 Rg on the equatorial plane 
and is elongated obliquely. The rarefied hot funnel region is 
shown by the red zone along the z-axis. 

On the other hand, when a is as large as 0.1, the shock 
wave formed initially in the inner region that expands out¬ 
ward and finally disappears from the outer boundary. At the 
later phases, a large scale convective motion is prevailed in 
the entire accretion flow. As a result, most of the input ac¬ 
creting matter is ejected from the outer boundary and the 
mass-outflow rate is modulated by a factor of few. 


3.2 The low angular momentum accretion flow 

The low angular momentum accretion flo ws without viscos¬ 
ity were examined in the previous works (lOkuda fe MoltenH 
I2012l : i0kudal2014i') . Using same sets of the model parameters 
of the specific angular momentum A and the specific total en¬ 
ergy e which were estimated from t he analysis of stellar wind 
of ne arby stars around Sgr A* (jCzern v fe Moscibrodzkal 
120081 : iMoscibrodzka, Das fc Czerml 20061 1, we showed that 
the irregularly oscillating shocks are formed in the inner re¬ 


3.3 The advection-dominated accretion flow 

Although the initial flow variables of case C are obtained 
by the self-similar solution of the advection-dominated flow, 
the evolution of the flow at the initial state is similar to 
that of case B. The upstream accreting gas falls towards 
the gravitational center and the outflows are generated from 
the surrounding region of the inner edge. These outflowing 
gas interacts with the downward accreting gas, resulting in 
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Figure 2. Subsequent snapshots of temperature contours at phases of t/10 6 = 2.899, 2.929, 2.960, 2.990, 3.005 and 3.036 (s) during an 
eruptive mass ejection in the curve of M 0 ut with a = 0.0 of Fig. f. The phases of the snapshots (a) and (f) are shown by two crosses in 
the solid line M ou t of Fig. 1 




Figure 3. Time variations of total luminosity L-y (erg s 1 ), mass- 
outflow rate Mout (gm s~ 1 ), mass-inflow rate M; n at the inner 
edge and shock position Rs/R s on the equatorial plane for case 
B with a = 10 -3 . 


Figure 4. Contours of ion temperature Tj with velocity vectors 
at t = 1.83 X 10® s of case B with a = 10 -3 . The evolutionary 
phase is denoted by “+’ in the curve M 0u t of Fig. 3. The velocity 
of light is indicated by the upper arrow. 


the formation of oblique shock. However, since the flow has 
larger angular momentum than that in case B, the shock on 
the equatorial plane continues to move outward gradually 
and finally reaches the outer disc boundary. As a result, the 
phenomenon of shock formation is never recovered. After a 
long integration time more than ~ 100 orbital period at the 
outer boundary, the entire flow seems to settle down in a 
steady state, although the luminosity is still modulated by 
a very small factor. On the other hand, the mass-outflow 
rate remains as large as 10% of the input accretion rate and 
varies by a factor of ten or so. 

Fig. 5 shows the time variations of the total luminosity 


Z/t erg s' 1 and the mass-outflow rate M out g s” 1 for case C 
with a = 0.1. In the plot, dashed horizontal line denotes the 
input accretion rate ~ 10 20 g s _1 . It should be noticed here 
that the mass-outflow ejected from the outer boundary sur¬ 
face with the polar angle £ ^ 60° is presented in Fig. 5 due 
to the following reason. We imposed free-floating conditions 
on the outer atmospheric boundary which allow only for out¬ 
flowing matter. However, recent HD simulations of the hot 
accretion flow on large scales , namely, from the black hole 
horizon to 10 Bondi radius (iLi, Ostriker fe Sunvaevl|2013l) 
indicate that the above boundary conditions may be invalid. 
This is because convective motions prevail even in the outer 
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Figure 5. Time variations of total luminosity Lx (erg s _1 ) and 
mass-outflow rate M ou t in case C with a = 0.1, where the dashed 
horizontal line denotes the input accretion rate. 



r/Rg 


Figure 6. Contours of density p(g cm -3 ) and velocity vectors at 
t = 5.0 X 10® (s) for case C. The velocity of light is indicated by 
the upper arrow. 

boundary region used in our simulations and accordingly 
there may exist outer flows with not only positive velocities 
but also negative ones at the outer boundary. On contrary, 
the boundary conditions adopted in this work seems to be 
justified in the polar region with £ ^ 60°, where the high 
velocity jets with only positive radial velocities appear in 
this conical region. 

In Fig. 6, we present the contours of the density at t = 
5.0 x 10 6 s for case C with a = 0.1. The corresponding 
contours of ion temperature are depicted in Fig. 7, where the 
velocity vectors are given by unit vectors. Here we find that 
the flow is unstable against convection and many convective 
cells are found in the inner region. It is remarkable that a 
torus disc with an inner edge of r ~ 3 R s and a concentric 
center at r ~ 6 R g is formed in the advection-dominated 
accretion flow. 

In Fig. 8, we show the variations of density p, ion tem¬ 
perature Ti, electron temperature T e and radial velocity v on 
the equatorial plane as function of radial coordinate. Note 



Figure 7. Same as Fig. 6 but for the ion temperature T; with 
unit velocity vectors, where many convective cells are found in 
the inner region. 


that the density and the temperatures on the equator have 
their maximum values at r ~ 6 R g and decrease sharply in¬ 
ward. At r ~ 3-Rg, the temperatures start increasing discon- 
tinuously although density continues to decrease even fur¬ 
ther. In Fig. 9, we zoom the inner edge of the torus disc to 
represent the contours of density with velocity vectors at the 
same phase as Figs. 6 and 7 for clarity. It is clear from the 
figures that the torus disc with the center at r ~ 6 R g has 
the inner edge of the disc at r ~ 3R g . We observe that the 
outflow begins at r ~ 3 R g on the equator and goes up along 
the outer surface of the torus disc. In the region r < 3R g , the 
gas falls toward the black hole and the density at the torus 
edge becomes very small and the flow has a discontinuous 
surface of the temperature there. 

In order to understand the irregular oscillation of mass- 
outflow, we examine the time-dependent animations of den¬ 
sity and temperature contours with velocity vectors. Then, 
we notice the formation of a series of hot blobs near the inner 
edge of the torus and their subsequent evolution along the 
outer surface of the torus. These hot blobs are originated in 
the innermost region between the inner edge ( r/R g ~ 3) of 
the torus disc and the inner edge ( r/R g = 1.5) of the flow. 

We examine another case with a = 10 -3 for the 
advection-dominated accretion flow. In this case, we again 
find the similar results as seen for a = 0.1, i.e., the large 
modulation of the mass-outflow rate, the small modulation 
of the luminosity and the torus disc formation. The mass- 
outflow rate varies irregularly by more than a factor of ten 
and roughly occurs in every several hours. A series of hot 
blobs is also generated in the innermost region and grows 
up along the boundary zone between the outflow and the 
outer surface of the disc. The torus disc has the center at 
r ~ 10-Rg and the inner edge of the disc at r ~ 6R g . The 
total luminosity is obtained as ~ 10 37 erg s _1 which is one 
order of magnitude larger than the case with a = 0.1. This is 
due to the fact that the input density at the outer boundary 
is much larger as seen from equation (11). 
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Figure 8. Variations of density p (g cm -3 ), ion temperature 7j, 
electron temperature T e and radial velocity v on the equatorial 
plane as observed in case C. 



r/Rg 


Figure 9. Same as Fig. 6, but near the inner edge, where at 
r ~ 3f?g the mass-outflow begins upward along the disc surface. 


4 KELVIN-HELMHOLTZ INSTABILITY AS 
THE ORIGIN OF THE UNSTABLE NATURE 
OF THE OUTFLOWS 

We have examined that the mass-outflows in the geometrical 
thick accretion flows such as cases A, B and C are generally 
unstable and the mass-outflow rates vary by more than a 
factor of few up to ten. In cases of A and B, the apparent 
origin of the outflow activity is ascribed to the formation 
of hot blobs and its evolution toward the outer boundary. 
We have examined such modulations of the mass-outflow 
rate in other simulations of 2D rotating accretion flows 
around a stellar-mass (IOMq) and a supermassive (10® Mq) 
black holes over a wid e range of input accretion rate s as 
10 -7 ^ M/Me ^ 10 -4 dOkuda, Teresi fe Molteiill2009l ). In 
addition, examining 2D simulations for the disc of SS 433 at 
super-Eddington luminosities, we found remarkable modu¬ 
lations of the accretion rate near the inner edge that results 
in the formation of recurrent hot blobs with high temper¬ 
atures and low densities at the disc plane. Based on this 
observation, we proposed that this may explain the massive 


jet e jection as well as the QPOs pheno mena observed in SS 
433 ilOkuda. Lipunova &: Moltenill2009l ). 

What is the intrinsic origin of the unstable nature of the 
outflows which is relevant to the hot blobs? In this respect, 
we consider shear instability or Kelvin-Helmholtz instabil¬ 
ity that may lead to the hot blobs and the winding waves 
of density. The shear instability and the Kelvin-Helmholtz 
instability can generally occur when there exists a velocity 
shear in a single continuous fluid or if there is a velocity 
difference across the interface between two fluids. Finally, 
the theory predicts the onset of instability and transition to 
the turbulent fl ow in fluids of differe nt densities moving at 
various speeds (IChandrasekhailll96ll ) . Certainly, in cases A 
and B, there exists an interface between the outgoing flow 
and the downward accreting gas as the shock surface with 
discontinuous velocities and densities and the hot blobs are 
formed along this interface and are evolved toward the up¬ 
per shocked region. For case C, the hot blobs exist along the 
interface between the outflow and the outer surface of the 
torus disc. Therefore, the Kelvin-Helmholtz instability gen¬ 
erated at the interface may be developed into the hot blobs 
and the winding waves. 

For the sake of simplicity, let us consider an extreme 
type of stratification, namely a plane parallel two-layer sys¬ 
tem in which a layer of fluid with velocity vi and density 
pi floats over another layer with velocity V 2 and density p 2 
where a gravitational force g acts vertically to the layers. 
Then, gravity waves propagate through the interface sepa¬ 
rating the two layers and grow in time. This leads to over¬ 
turning in the vicinity of the interface and mixing over a 
height SH along the z-axis vertical to the layers, which is 
called the Kelvin-Helmholtz instability. In this simple pic¬ 
ture, the growth of the Kelvin-Helmholtz ins tab ili ty is deter¬ 
mined by a crit erion of Richardson number llChandrasekhail 

Il96ll: Ishul[l99^) . 


= -g{dp/dz) 1 
“ p{dv/dz ) 2 4' 


(14) 


Then, we have the Richardson number approximately as 


Ri — 


-g8H | p 2 - pi 
p(v 2 — Vl) 2 


(15) 


In Fig. 2, we find the oblique shock front present in 
the accretion flow. Regarding the shock front and 8H as 
the interface mentioned above and the shock thickness, re¬ 
spectively, we estimate the Richardson number. Across the 
interface, if the velocity difference parallel to the shock sur¬ 
face is not large enough, the criterion described in Eq. (14) 
is never satisfied. Fig. 10 illustrates the shock surface (thick 
line) as in Fig. 2 and the geometry of the shock front with 
thickness 8H where S(r,() and T(r, ( + 8() are the inter¬ 
section points of the shock front with a circle of radius r 
centered at the origin, ui is the angle between the gravita¬ 
tional force direction and the tangent to the shock surface 
at point S, r/ is the angle between the tangent at point S 
and x-axis. With this, the Richardson number in units of 
the Schwarzschild radius and the speed of light is given by, 

D _ 1 \P2-P1 I^Csinwcos (w-f) 

'0 — 7T z 1 To I (lb) 

2 p r ( V 2 — vi) 2 


where v\ and V 2 are the components of pre-shock and post¬ 
shock velocities parallel to the tangent at S. The physical 
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Figure 10. Geometry of shock surface with shock front thickness 
SH found in Fig. 2, where S(r, () and T(r, C + <5C) are the intersec¬ 
tion points of the shock front with a circle of radius r centered at 
the origin, u) is the angle between the gravitational force direction 
and the tangent to the shock surface at point S, 77 is the angle 
between the tangent at the point S and x-axis. 


Figure 11. Plot of density p (g cm -3 ), radial velocity v r and 
azimuthal velocity v v as function of polar angle ( at r/R g = 67.4 
on the shock front in Fig 2a, where the shock transition is found 
in the region of 63.9° ^ C ^ 67.5° and ~ 3.6°. 


variables at the points S and T are determined approxi¬ 
mately from the numerical data of the flow. As an exam¬ 
ple, we take a point on the shock front in Fig. 2a., whose 
radius is 67.4_R g . For this point, in Fig. 11, we show the 
variation of density p, radial velocity v T and azimuthal ve¬ 
locity v v as function of polar angle £ measured from the 
equatorial plane. The shock transition is found here in the 
region of 63.9° Sj £ ^ 67.5° and 5C, ~ 3.6°. From the 
numerical data, furthermore, we know pi, p 2 , (v r )i, (v r ) 2 , 
(tV)ij ( v v>) 2 , w(= r) - C), vi = (u r )icos 0 J + (u c )isin w, 
V 2 = (u r ) 2 COS (w — <5() + (v£ ) 2 sin (co — and finally estimate 
Ri. For case C, we do not find the shock interface found in 
cases A and B. However, there exists a discontinuous bound¬ 
ary zone between the torus disc and the outflowing gas where 
change of temperature and velocities occur very sharply as 
observed in Figs. 8 and 9. We apply the above method at 
this boundary zone to obtain the Richardson number. 

In Fig. 12, we plot the distribution of the Richardson 
number Ri as a function of radial distance r/R g obtained for 
the oblique shocks in Fig. 2(a) of case A (open squares) and 
in Fig. 4 of case B (diamonds), and for the boundary zone 
around the torus disc as in Fig. 9 of case C (closed circles). 
This shows that the criterion (Eq. 14) for Kelvin-Helmholtz 
instability is not satisfied near the equatorial plane but is es¬ 
tablished apart from the equator for all cases. This suggests 
that the unstable property of the outflows would possibly 
be ascribed to the Kelvin-Helmholtz instability. 


5 A POSSIBLE MODEL TO THE FLARES OF 
SGR A* 

The supermassive black hole candidate Sgr A* possesses 
a quiescent state and a flare state. The flares of Sgr A* 
have been detected in multiple wavebands from radio, sub¬ 
millimetre, IR to X-ray. The unstable mass-outflow activi¬ 
ties examined in this work demonstrate both small and large 
amplitude modulations of the total luminosity and the mass- 



Figure 12. Variation of the Richardson number R; as function of 
radius on the shock surface in Fig. 2(a) of case A (open squares), 
in Fig. 4 of case B (diamonds) and on the boundary zone around 
the torus disc in Fig. 9 of case C (circles). The corresponding fit¬ 
ted curves are shown by thick solid, dash-dotted and long-dashed 
lines, respectively. The horizontal dashed line denotes the critical 
value Rf r = 0.25 for the Kelvin-Helmholtz instability. 


outflow rate which seems to be relevant for explaining the 
flaring phenomena observed in Sgr A*. The total luminosity 
found in cases B and C satisfactorily explains the existence 
of low luminosity of ~ 10 36 erg s -1 for Sgr A*. Similarly, 
the irregularly eruptive mass-outflows may correspond to 
the flare activities of Sgr A*. 

The observations of Sgr A* show that the duration of 
the X-ray and IR flares are typically of 1-3 hours and the 
flare events usually occurs few times per day. In the sim¬ 
ulation of cases B and C, we find that the duration and 
the interval between the successive eruptive mass-outflows 
roughly are in agreement with the flares of Sgr A*. But, 
the small modulations of luminosity does not correspond 
to the large amplitudes of the flare intensities observed 
at various wavelengths, because the amplitudes at radio, 
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IR and X-ray are ro ughly varied by f a ctors of 1/2, 1-5 


and 45, respectivelv (Genzel et al., 2003 

: Ghez et alJ 

2004 

Eckart et al. 20061: Mever et al. 2006al.b: 

TriDDe et al. 

2007; 

Yusef-Zadeh et al. 2009. 2011|). It is well known that the sub- 


milimeter flare lags the X-ray flare roughly by two hours. 
From these facts, we propose a possible model of the unsta¬ 
ble mass-outflow applicable to the flare phenomena of Sgr 
A*. This is basically a scenario of the temporal sequence of 
disc-jet coupling observed in the X-rays, IRs a nd radio wave¬ 
lengt hs of the microquasar GRS 1915+105 dMirabel et al.l 
1 199811 . Here, we consider that the total emission in our mod¬ 
els corresponds to the permanent emission of the quiescent 
state of Sgr A*, whereas only the emission of Sgr A* flare 
is produced at the distant region far away from the Sgr A* 
through the interaction between the surrounding dense in¬ 
terstellar matter and the high-velocity jets. These jets are 
originated in the unstable mass-outflows as examined here. 
The ejected wind carries sufficient outflowing mass of rate 
~ 10 1 ® - 10 20 g s -1 and a total kinematic energy of ~ 10 36 
-10 37 erg s _1 to reproduce the flare phenomena. However, 
one needs to incorporated other mechanisms also so that the 
high-velocity wind could be accelerated up to the relativis¬ 
tic speed. This is because the outflow velocities at the outer 
boundary in cases B and C are limited to sub-relativistic 
regime as ~ 0.01 — 0.1c. 

When we consider the flare activities of Sgr A*, it is 
required for us to treat MHD model because it is widely be¬ 
lieved that the existence of the magnetic field is inevitable 
from the spectral analyses of Sgr A* and that magnetic tur¬ 
bulence induced by magneto-rotationa 1 instability is respon¬ 
sible for the angular momentum transfer. In this respect, by 
analogy with the coronal mass ejection on the Sun, Youan 
et al. (2009) proposed a MHD model for the formation of 
episodic jets from black holes associated with the closed 
magnetic fields in an accretion flow. In their model, shear 
and turbulence of the accretion flow deform the field and re¬ 
sult in th formation of a flux rope in the corona and the flux 
rope is ejected through the magnetic reconnection process 
in a current sheet under the magnetic field. 

Although the unstable outflows examined here are re¬ 
stricted to the non-magnetized accretion flow, we notice 
that the comparative studies of 2D hydrodynamical and 
magneto-hydrodynamical accretion flows from a same thick 
torus around black holes show many similarities such as 
highly time-dependent accretion rate and Kelvin-Helinholtz 
instability in spite of t he different mechanism of angu¬ 
lar momentum trans fer (iStone, Pringle fc Begelmanl 1 19991 ; 
IStone fc Pringlell200ll l. They also show that, for strong mag¬ 
netic field, a small fraction of outflows near the pole is es¬ 
caped as a powerful MHD wind, whereas the outflows in 
the hydrodynamical simulations are bound. This suggests 
that the unstable nature of the outflows from the thick ac¬ 
cretion flow may be furthermore enhanced in 2D magneto- 
hydrodynamical simulation. 


6 SUMMARY AND DISCUSSION 

We examined the mass-outflow in the hot geometrically 
thick accretion flows such as the low angular momentum 
accretion flow and the advection-dominated accretion flow 


around black holes, using 2D time-dependent HD calcula¬ 
tions. The results are summarized as follows. 

(1) In the low angular momentum flows, the inward 
accreting matter on the equatorial plane interacts with the 
outward moving gas and the centrifugally supported oblique 
shock is formed along the interface of both the flows provide 
the viscosity parameter a is as small as a ^ 10~ 3 . The hot 
blobs, which lead to the eruptive mass-outflow, are gener¬ 
ated at the inner shocked region and grow up toward the 
outer boundary. 

(2) The advection-dominated accretion flows attain fi¬ 
nally to a torus disc with the inner edge of the disc at 
3 ^ r/R g ^ 6 and the center at 6 R g ^ r Sj 10R g . The 
hot blobs are formed near the inner edge of the torus disc 
and grow up along the outer surface of the disc. 

(3) In both models, the mass-outflow is unstable and 
eruptive in nature. As a result, the luminosity and the mass- 
outflow rate are modulated irregularly and especially the 
relative modulation amplitudes SM ou t/M ou t of the mass- 
outflow rate are as large as a few to ten. 

(4) The origin of the unstable mass-outflow is ascribed 
to the Kelvin-Helinholtz instability which is ascertained by 
examining the Richardson number criterion (Eq. 14) for 
Kelvin-Helinholtz instability at the inner region of the flow. 

(5) The flare phenomena of Sgr A* may be explained 
by the unstable mass-outflow examined here, together with 
the scenario of the temporal sequence of disc-jet coupling 
observed in the X-rays, IRs and radio waveleng ths of the 
microquasar GRS 1915+105 dMirabel et al.lll998f ). 

The unstable mass-outflow seems to be inevitable as far 
as we consider the thick accretion flows with small angular 
momentum around the black holes. The present study shows 
that both the models of the low angular momentum accre¬ 
tion flow and the advection-dominated accretion flow pro¬ 
duce the eruptive mass-outflows. These models would repro¬ 
duce the reasonable spectra fitted to the observed ones of Sgr 
A* if we take in the appropriate model parameters, such as 
the accretion rate M. the viscous parameter a, the magnetic 
field strength /3 and a hybrid electron population of thermal 
and nonther mal particles, fo l lowing the i dea and m ethod 
described bv lYuan, Quataert fe Naravanl d2003l . I20041) . Fur¬ 
ther, the evidence of angular momentum would confirm the 
validity of the models. Moreover, it is required to perform 
magneto-hydrodynamical research instead of only the hy¬ 
drodynamical one. This is because the magnetic field will 
play important roles not only in the small-scale region of 
the accretion flow but also the large-scale of the interstellar 
space. 
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